Damping of Growth Oscillations 
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Computer simulations and scaling theory are used to investigate the damping of oscillations during 
epitaxial growth on high-symmetry surfaces. The crossover from smooth to rough growth takes 
place after the deposition of (D/F) s monolayers, where D and F are the surface diffusion constant 
and the deposition rate, respectively, and the exponent S = 2/3 on a two-dimensional surface. 
At the transition, layer-by-layer growth becomes desynchronized on distances larger than a layer 
coherence length proportional to I 2 , where I is a typical distance between two-dimensional islands 
in the submonolayer region of growth. 

PACS numbers: 81.10.Aj, 81.15-z, 68.55-a, 05.70.Ln 



I. 1 INTRODUCTION 

Layer-by-layer or Frank-van der Merwe growth |l| is 
a growth mode observed in molecular beam epitaxy and 
other deposition methods which allows precise control of 
chemical composition of layers down to atomic thickness. 
It is therefore particularly well suited for the fabrication 
of novel electronic devices. 

The key microscopic processes in layer-by-layer growth 
are deposition of atoms onto a high-symmetry surface 
and diffusion of adatoms on the surface. The adatoms 
meet and form dimers which then grow into islands of 
monoatomic height whose edges capture most of the 
adatoms during the deposition of one monolayer. When 
the island edges become less available due to coalescence 
of islands, formation of dimers and islands in the next 
layer begins. The density of atomic steps - and all other 
quantities sensitive to the surface morphology - thus os- 
cillates in time. 

Generically, these oscillations are damped: Layer-by- 
layer growth is only a transient. Possible reasons Q in- 
clude (i) cessation of periodic formation of islands on the 
surface and transition to step flow growth or (ii) 

roughening of the surface ||. If the substrate tempera- 
ture is increased, damping becomes stronger in the first 
and weaker in the second case, allowing to discriminate 
between the two. Ignoring the possibility of inhomoge- 
neous deposition (cf. Ref. JjJ), surface roughening can 
have two different sources. If intcrlayer transport is in- 
hibited by step-edge barriers, one obtains the growth in- 
stability predicted by Villain ||. If no such instability 
occurs, the surface may still roughen due to fluctuations 
in the intensity of the deposition rate. Only the latter 
case is considered in this paper. 

The question of how long layer-by-layer growth persists 



is of immediate practical importance. If the answer is 
known, one can devise and optimize annealing schedules 
for growing thicker films while maintaining a smooth sur- 
face. For this purpose, it is important to know both the 
damping time t and the length scale over which layer-by- 
layer growth is synchronized. This layer coherence length 
£ is a new characteristic length that determines, e.g., the 
annealing time needed to reestablish a flat surface || be- 
fore growth can be continued. 

A theory has been proposed recently |n],[llj which pre- 
dicts that the damping time and the layer coherence 
length depend on the typical distance £ between islands 
in the submonolayer region of growth as 

Ftoc£ M /^ and loc^^. (1) 

In this paper, we present the theory of oscillations 
damping and detailed numerical evidence of the valid- 
ity of its predictions based on extensive computer sim- 
ulations of a minimal, one-parameter model at surface 
dimension d = 2. The methods of extracting the damp- 
ing time and the layer coherence length from the surface 
morphology evolution are outlined and thoroughly dis- 
cussed. 

The characteristic length £ (and thus also t and £) has 
a power-law dependence on the ratio D/F of the surface 
diffusion constant to the deposition rate: 

£oc(D/Fy (2) 

(see and references given in |^||). The exponent 7 
depends on the dimensionality d of the surface and the 
(possibly non-integer) dimension dj of the islands. It also 
depends on whether or not desorption of adatoms or dif- 
fusion of dimers or larger clusters is negligible. Finally, 
7 is a function of the critical cluster size i* for the for- 
mation of a stable nucleus. For the case considered here 
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in more detail (d = 2, df = 2, i* = 1, no desorption, 
immobile clusters), the value is 7 = 1/6. 

The layer coherence length £ as well as the damping 
time t play the roles of natural cutoffs in the continuum 
growth equation at small length and time scales. For 
t > t one expects that the surface exhibits self-affine 
scaling [Ol: 



«;(*)- OL(£(t)/I) 



C 



and 



ffi - t (t/t 



(3) 



Here w is the root mean square variation of the film thick- 
ness (the surface width) , a± the thickness of one atomic 
layer (which, for convenience, is set to one in the fol- 
lowing), and £; the correlation length up to which the 
surface roughness has fully developed at time t. £ is the 
roughness exponent and z the dynamical exponent. The 
dependence of £ and t on the microscopic growth param- 
eters will be derived next. 



II. 2 THEORETICAL RESULTS 

Coarse-graining the surface configuration at a given 
time over a length scale of the order of £ , one can write 
down an evolution equation for the variable h(x, t). Since 
particle desorption can be neglected under conditions 
typical for molecular beam epitaxy, the equation can be 
written in the form of a conservation law, 



d t h(x,t) — —Vj(x,t) +r](x,t). 



(4) 



j is the surface diffusion current, and r\ is white noise 
with second moment 



(r)(x,t)r)(y,s))=F6 d (x-y)6(t-s), 



(5) 



which describes the fluctuations in the deposition rate. It 
was proposed by Villain ||] that in growth processes far 
from equilibrium where local chemical potentials along 
the surface are ill defined, diffusion currents should be 
driven by gradients in the growth-induced, nonequilib- 
rium adatom density n 1l4| , 



(6) 



On a singular surface, the balance between deposition 
and capture of adatoms at steps leads to a stationary 
adatom density n = n a of the order of |ll| n — (F/D)£ 2 . 
On a vicinal surface, the adatom density is reduced due 
to the presence of additional steps. However, this effect 
is felt only if the miscut m — \ Vh\ exceeds l/£, in which 
case n ~ (F/D)m~ 2 . In terms of a coarse-grained de- 
scription of the surface this implies that the local adatom 
density depends on the local miscut or surface tilt. A 
useful interpolation formula which connects the regimes 
m < l/£ and m > l/£ is pi 



i(V/i) = 



n 



1 + (£\7h) 2 

(F/D)£ 2 - (F/D)£ 4 (Vh) 2 



(7) 
(8) 



Inserting the leading quadratic term of this gradient ex- 
pansion into (^), which is appropriate for describing long- 
wavelength fluctuations around the singular orientation, 
one obtains 



with 



VA(V/i) 2 



A = Ft 



(9) 



(10) 



Considering Eq. (jij) and Eq. (JsJ) one sees that the phys- 
ical dimension of A is (length) F/ (time-height). Within 
the continuum description, the only characteristic length 
and time scales are the layer coherence length and the 
damping time, whereas the lattice constant aj_ has been 
chosen as a unit of height. Therefore 



A oc £ A /t 



(11) 



on dimensional grounds. 

Finally, the number of particles deposited during the 
time t onto an area £ d is Ft £ d ± {Ft £ d ) Y l 2 . Thus the 
fluctuatio n of t he film thickness over the distance £ is 
w(t) ~ VFt £ d /£ d . At t this should be the thickness of 
about one atomic layer, w(t) ~ 1, which results in 



Ft = 



(12) 



Combining fll0|), fllj) and (|12J) one obtains Eq. fly), or 
Ft~{D/F) & and £~{D/F) s,d (13) 
with the exponent 



= 7. 



(14) 



Notice that the layer coherence length £ is substantially 
larger than the characteristic distance £ between islands 
{£~£ 2 atd = 2). 



III. 3 MODEL 

In our model, atoms are deposited onto the (100) sur- 
face of a simple cubic lattice with the rate of F atoms per 
unit time and area. The surface size is L x L = 128 x 128 
Jl7j . Atoms with no lateral neighbors are allowed to dif- 
fuse with diffusion constant D. Atoms with lateral neigh- 
bors are assumed to be immobile so that, e.g., dimers are 
immobile and stable. Growth commences on a flat sub- 
strate, h(x, 0) = for all sites x. On deposition at x, 
h(x, t) is increased by one. We neglect barriers to inter- 
layer transport (Ehrlich-Schwoebel barriers H| ) so that 
the only parameter of the model is the ratio D/F. 
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FIG. 1. Semilogarithmic plot of the upper envelopes of 
the oscillations of the kinematic intensity / as a function of 
time. The minima of the oscillations are zero, as indicated 
schematically for D/F — 500. Five different D/F values are 
evaluated, spanning almost two decades. The dotted line de- 
notes the value I — 0.05 used for the determination of the 
damping time. 



IV. 4 DAMPING TIME 

First we present the results for the damping time ex- 
tracted from kinematic intensity data, 



/ = ((N cvcn - N odd ) 2 /L 2 ) 



(15) 



(see Fig. [|; N cvcn (N 0< id) denotes the number of atoms in 
even (odd) layers). The brackets (. . .) denote averaging- 
over different runs. The same analysis was done for the 
surface width with equivalent results for the damping 
time. 

The kinematic intensity oscillates between zero and 
maxima which decrease until the oscillations vanish. We 
measure t as the time where the maxima of the kinematic 
intensity drop below / = 0.05 p9| . The results are shown 
in Fig. |. 

Obviously there are strong corrections to scaling which 
can be attributed to an offset to > 0: 



(16) 



to plays the role of a cutoff for the validity of our scal- 
ing theory. (D/F)q = (to/Af) 1 / 6 can be interpreted as 
that value of D/F below which the oscillations are not 
observable anymore pp| . 

A three-parameter fit to the data shown in Fig. || gives 
an exponent 



6 = 0.69 ±0.05, 



(17) 



where the error bar reflects the variations obtained when 
the evaluation method is modified jl9],^0| or when the 
data for the surface width are evaluated in the same way. 
This value is in good agreement with the theoretical pre- 
diction of 5 = 2/3 (see Eq. Q) for compact islands. 
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FIG. 2. Damping time t determined from Fig. as a func- 
tion of D/F. The straight line has a slope 2/3, in agreement 
with the theoretical prediction. The dotted line is the best fit 
according to Eq. 
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FIG. 3. Scaling the deposition time with the damping 
time [Eq. leads to a good data collapse of the curves 

presented in Fig. [l]. 

In our simulation, neither detachment of adatoms from 
islands nor edge diffusion are considered. Therefore, the 
islands are fractal for large diffusion lengths ]2l[-p3]| with 
the fractal dimension df ~ 1.72 of two-dimensional dif- 
fusion limited aggregation p4| . Then 7 changes from 
7 = 1/6 to 7 ~ 0.175. This leads to theoretical pre- 
diction S ~ 0.70 which is also within the error bars of 
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Eq. (J!?]) p5|| . However, in the present case the values of 
D/F are sufficiently small, so that this complication may 
be ignored. 

The scaling plot of the kinematic intensity with the 
time divided by the damping time according to Eq. ( |l6| ) 
(see Fig. ||) confirms the validity of the approach used. 



V. 5 LAYER COHERENCE LENGTH 

The measurement of the height difference correlation 
function [t25| 



A\t) = (w 2 (t)) %t+T] - Mi))f M+r] 



(19) 



G{X, t) EE ([h(x , t) - h{x + X, t)} 2 }, 



(18) 



evaluated at t = i shows that G has a maximum. This 
can be explained as follows. At t — t the probability to 
find the surface at the same height as at a reference point 
xq is minimal at a distance x — xq corresponding to the 
layer coherence length. For larger distances, deviations 
from the average height are essentially uncorrelated. At 
very small distances, their correlation is positive, while 
around £ they are anticorrelated. This is how the data 
denoted by the squares in Fig. || were obtained. The 
result is in good agreement with the predicted exponent, 
cf. Eq. (13). (Note that this method could be also used 
for experimental determination of £.) 
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FIG. 4. £ as a function of D/F, measured using the height 
difference correlation function (squares) , and finite-size analy- 
sis (circles) , respectively. The solid lines have the theoretically 
predicted slope 1/3 for compact islands, see Eq. (|l3|). 

An alternative method of measuring £ is to carry out a 
finite-size analysis in the following way. The surface does 
not roughen when the linear system size L is smaller than 
£. In this case, the amplitude of the growth oscillations 
becomes stationary after a transient time, and the oscil- 
lations never die out. We monitored the variance 



of the surface width during the layer completion time 
t ee 1/F, where (. . .}u : t+r] denotes the time average over 
the interval [t, t + t}. Its stationary value decreases with 
increasing system size and ultimately becomes equal to 
the statistical fluctuations of w when the system size is 
big enough so that the oscillations can die out completely. 
The values of £, denoted by the circles in Fig. [l] represent 
the linear system size L at which the stationary value 
of A(t) drops below 0.37. Both methods of measuring £ 
are in excellent agreement with each other and with the 
theoretical prediction pa]. 



VI. 6 CONCLUSIONS AND OUTLOOK 

We have presented a theory for the damping of growth 
oscillations caused by kinetic roughening. We have shown 
that the results of numerical simulations of a minimal 
model compare very favorably with the theory, and di- 
rectly determined two key quantities, the damping time 
and the layer coherence length. The instability associ- 
ated with barriers to interlayer transport ]^| may com- 
pete with the kinetic roughening mechanism as a source 
of oscillation damping |2(|. This, as well as the transi- 
tion to step-flow growth on vicinal surfaces will lead to 
different power laws for the damping time and the layer 
coherence length. This remains for future research. 

The results of this paper can be directly verified by 
diffraction or real-space surface sensitive techniques. The 
determination of the damping time and, in particular, of 
the layer coherence length as a function of growth con- 
ditions should be possible using the methods outlined 
above. 
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